Orientation field model for chiral branching growth of bacterial colonies 
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We present a new reaction-diffusion model for chiral branching growth of colonies of the bacteria 
Paenibacillus dendntiformis. In our model the bacteria are represented by a density field with non- 
linear diffusion and a complex scalar field which represents bacterial orientation. The orientation 
field introduces anisotropy into the fiux of bacteria, representing self-propulsion along their long 
axes. The model can also reproduce tip-splitting growth of other strains (shorter bacteria) of the 
same species. The model can capture changes of small number of bacteria, thus it can be used to 
study the open question of transitions between tip-splitting and chiral dendritic growth. 



Bacteria display various chiral properties. Mendelson 
et al. showed that long cells of Bacillus subtilis can 
grow in helices, in which the cells form long strings that 
twist around each other. They have also showed that the 
chiral characteristics affect the structure of the colony. 
Matsuyama et al. ^ have found that colonies of B. sub- 
tilis can grow into tip-splitting patterns with a global 
rotation about the center of the colonies (global twist). 
Ben- Jacob et al. have found that some strains of 
Paenibacillus dendritiformis (see Ref. [|] for identifica- 
tion of the bacteria) exhibit similar patterns with global 
twist. They have also found that other strains (referred 
to as C (chiral) morphotype) present different chiral prop- 
erty, chiral dendritic growth with local twist (which they 
refer to as strong chirality). 

Colonial patterns of C bacteria grown on semi-solid 
agar are characterized by chiral dendritic branching pat- 
terns, where the branches are narrow and are twisted 
with the same handedness (Fig. |l]). All colonies of this 
bacteria grown in similar conditions show the same hand- 
edness. Side branches are usually emitted to the convex 
side of the arced branches (see Refs. for morpholo- 

gies and studies of C bacteria) . 




FIG. 1. Chiral colonial growth of Paenibactllus dendri- 
tiformis, strain C. a) Global view of a colony shows thin 
branches, all twisted with the same handedness. Colony is 
grown at 2g/l peptone level and 1.25% agar concentration, 
b) Optical microscope observations of branches of C colony. 
x20 magnification of a colony at l.&g/l peptone level and 
0.75% agar concentration, the anti-clockwise twist of the thin 
branches is apparent. The curvature of the branches is the 
same throughout the growth. 

2D chiral branching patterns are also observed in non- 



living systems. Shapes resembling Sea Horse, or S, 
are formed during deposition of thin films of fullerene- 
tetracyanoquinodimethane (C60-TCNQ) or pure TCNQ 
[Q. Modeling of the system indicates that the apparent 
curvature of the branches is actually a strong bias in se- 
lection of splitting branches, and the branches themselves 
are not curved. The patterns of bacterial colonies share 
more resemblance with patterns formed during compres- 
sion of monolayers of various chiral molecules at the air- 
water interface |^ , and electro-chemical deposition under 
a uniform magnetic field |9| . Various modeling attempts 
indicates that the processes forming these patterns are 
related to processes of solidification (see Refs. [ p^p] ] for 
the differences between the processes of solidification vs. 
colonial growth). A different model is required for chiral 
branching growth of bacterial colonies. 

The Communicating Spinors model of Ben- Jacob et al. 
is an atomistic (discrete entities) model which de- 
scribes chiral branching growth of bacterial colonies. The 
model presented here is not a mean field model of the 
spinors model. The two models complement each other 
and highlight different biological features. Each spinor 
in the spinors model represents a large groups of about 
10^-10'* bacteria. This coarse graining makes simulations 
computationally feasible, but prevents modeling of pro- 
cesses which cannot be averaged over large groups (specif- 
ically, single-cell events). Continuous model is more ap- 
propriate for this purpose (when interpreting densities as 
probability densities) . The model presented here can 
be used to study the effect of mutations and transitions 
on colonial morphologies, as will be shown elsewhere. 

Detailed description of the materials and methods used 
in the experiments can be found in Ref. ||l|l. The bac- 
terial colonies are grown on top of agar (semi-solid jelly) 
with peptone as a nutrient. During colonial growth bac- 
teria are confined to within a layer of lubricating fluid, 
which is extracted from the agar by the bacteria them- 
selves. Inside the layer the bacteria propel themselves 
in a random-walk like motion. Optical microscope ob- 
servations indicate that the length of C bacteria is up to 
50 times their width. In narrow branches, bacteria are 
aligned with their neighbors, and their movement is quasi 
ID random walk along their long axis. 
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In order to capture the details of bacterial dynamics 
we define the bacterial density per angle b{x, 9, t), where 
X g IR^ is a position and G [0, tt] is an angle of orien- 
tation. We also use the bacterial density B{x,t), which 
is defined as the mean of b over 9. The equation for the 
bacterial dynamics is 

b =^ V ■ [Db{B)-Doi9)Vb] - deJe + Gib,n) (1) 

where V operates in the spatial dimensions, Di, is a non- 
constant diffusion coefficient , Jg is the flux in the angle 
dimension which represents changes in bacterial orienta- 
tion, G is a reaction term representing growth and death, 
and n is the nutrient concentration. We take the diffusion 

matrix Do (9) = R(6l)^ ^ R{9), where R(0) is 

a rotation matrix. D^i is a constant coefficient for diffu- 
sion in the bacterial direction 9 (due to self-propulsion 
forward and backward) and Di < Dd is a constant coef- 
ficient for diffusion in the lateral direction, due to fluctu- 
ations. Eq. (|^) is invariant under rotation in space and 
translation in 9 by the same angle. 

It was shown pj] , p^ -p^ that bacterial movement in 
a self-produced layer of fluid can be approximated by a 
non-linear diffusion, where the diffusion coefficient is pro- 
portional to the bacterial density to a power greater than 
one. The proportion constant is related to the agar dry- 
ness through the rate of absorption of the fluid into the 
agar. Hence we take Db{B) = [k > 2), where the 
proportion constant is included in and Di. We will 
assume that bacterial friction with the agar is propor- 
tional to the agar dryness, and as the bacterial velocity 
is inversely proportional to the friction, and Di are 
proportional to the agar dryness to the power —2. 

Following ||f^,|ll| we take a simple form for the growth 
function, G{b,n) = 6(n — /i), where > is a rate of 
conversion into immobile sporulating cells. A more accu- 
rate form should have saturation for high values of n/B, 
but the linear form is a reasonable approximation for low 
levels of nutrients, as is the case in the bacterial colonies. 

Bacteria change their orientation (Jg) in response to 
their neighbors orientation. It is convenient to use an 
auxiliary complex orientation field p, which is defined as 
p{x, t) — ^ Jq b{x, 9, i)e~'^^ d9. The mean orientation of 
the bacteria is along the vector (±Re^, ^luiy/p) . 

We found that it is sufficient to track the dynamic of 
B and p, instead of the full dynamic of b. Using Fourier 
expansion of Eq. (|l|) in 9 and taking the first two re- 
sulting equations we derive the Orientation Field (OF) 
model, which includes equations for B, p, the nutrient 
concentration n and the density of immobile (sporulat- 
ing) bacteria s (at this stage s is used only to record the 
history of the colonial development): 

B^V ■{b'' [DiVS + 2Re(D2Vp)]} +B(ra-/i) (2) 
P = V-{b'' [D;VB + DiVp]} + p{n - p) 



+a{B,\p\)p + ^{B''-^VB,p)ip (3) 
n = D,y'^n - nB (4) 
s== fiB (5) 

where n is scaled to units of bacterial mass and £)„ is its 
diffusion coefficient, a and 7 are real valued functions, a 
decomposition into orthogonal elements of the derivation 
of dgJe (we give here the functions a and 7, not the orig- 
inal Je). Here Di = and D2 = ^"'^' M , where 

M = ^ ^ (from here on we take all constants to be 

real, unless otherwise stated). Eqs. (||j|) are invariant 
under a rotation by (j) and a multiplication of p by e"'^*^. 
For the co-alignment function a we take 

a(B, \p\) = ~ADe + u [B - \p\)l{B + \p\) (6) 

The first term in the RHS results from linear diffusion of 
6 in the 9 dimension (with Dg the diffusion coefficient) 
and the second term is an alignment of non-align bacteria 
with the mean orientation (with v being the rate of this 
co-alignment). The exact form of a is not important, as 
long as a(|p|) has at most one positive root with negative 
derivative in the range [Q,B] and a(|p|) in non-positive 
outside the range (0,5). 

The rotation function 7 results from the only process 
that brakes left-right symmetry, the bacterial bias in their 
tumbling when they are placed in a thin layer of liquid. 
The bacteria are restrained by their neighbors, and the 
restraints on rotation are weakest near the boundary of 
the branches. Due to the non- linear diffusion ]ri| , p^ -p7[, 
the branch's boundary is identifiable by the divergence 
of VB (for the initial conditions specified below). The 
term B^~^V B is always finite and can be used to identify 
the boundary (as p is bound by -B, it can replace B in 
the coefficient). 7 can be written as a function of A, the 
angle between the branch's boundary and the bacterial 
mean orientation: 

7(a, (3) = Re [7° - 7"/? {a^Ma) / |a|] 

= \!3a\ [2Re(7") sin^ (A) -I- Im(7°) sin (2A)] (7) 

where 7'^ is a complex constant. Re (7*^) is a measure of 
the anti-clockwise bias of the bacteria at the tip of the 
branch and Im(7'^) is a measure of the torque aligning 
the bacteria in parallel to the boundary. The rotation at 
the tip of the branch is also restrained by friction with 
the agar, and from the same reasoning that related Dd 
and Di to the agar dryness, we deduce that Re (7*^) is 
inversely proportional to the agar dryness. 

As initial conditions (in a circular 2D geometry), we set 
n to have uniform distribution of level no, -B to have com- 
pact support at the center where it is positive, and the 
other fields to be zero everywhere. We solve the model 
numerically using a 2nd order explicit scheme. In order 
to reduce the implicit anisotropy of the scheme, we use 
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tridiagonal lattice and multiply the bacterial diffusion co- 
efficients by a quenched noise with mean 1. The quenched 
noise represent inhomogeneities of the agar surface. We 
show in Fig. ^ that the model can indeed reproduce the 
microscopic bacterial dynamics and the chiral branching 
patterns with the local twist. 





Results of numerical simulation of the OF model, 
a) Global view of a simulated colony with a local twist of 
branches. Densities of B + s are indicated by gray levels. Pa- 
rameters values are fc = 3, Dd ~ 0.0625, Di — 0, fj, = 0.1, 
De = 0.001, u = 0.5, 7° = 0.0075 - i0.006, L»„ = 1 and 
no = 1.5. b) A close look at the pattern of (a), showing the 
details of the curved branches. 



The Non-Linear Diffusion model (NLD) for colonial 
growth can reproduce tip-splitting branching patterns 
[ pd] , p^ -[l7[ of a related morphotype, T morphotype, with 
shorter bacterial cells. The OF model reduces to the 
NLD model if the self-propulsion is not primarily along 
the bacterial long axis {Dd — Di), if the bacteria do not 
tend to co-ahgn {f < iDg or {v - iDo) <€. {ly + ^Dg) 
), or if bacteria at the tip of the branch tend to rotate 
too freely (|7°| {v — ADq) > ID2I). In all these cases 
|Re(D2Vp)| < l-DiVBl everywhere and the OF model 
produce tip-splitting branching patterns. 

The response of the simulated growth to initial food 
concentration and agar dryness is shown in Fig. ||. 
In agreement with experimental observations the 
main effect of the parameters is on the global density 
of branches, while the curvature of the branches is only 
weakly affected. 






FIG. 3. The response of simulated colonies to experimen- 
tal control parameters. Figures (a) and (b) present change in 
initial food concentration with no = 1 and no = 2 in (a) and 
b) respectively. All other parameters are the same as in Fig. 

The expansion time of the colonies, normalized by the 
expansion time in Fig. |2|a, is 2.9 and 0.4 in (a) and (b) respec- 
tively. Figures (c) and (d) present change in agar dryness A. 
A is related to the model's parameters with Re(7'^) = 0.009/A 
and Dd = 0.09/ A^. In (c) and (d) A = 1 and A = 1.5 re- 
spectively, with all other parameters are the same as in Fig. 
pk, where A — 1.2. The normalized expansion time of the 
colonies is 0.8 and 1.2 in (c) and (d) respectively. 



It was shown |l8|,|,|l|,|l that many features of the colo- 
nial patterns are explained when food chemotaxis and 
chemotactic signaling (a chemotactic response to a ma- 
terial emitted by the bacteria themselves) are modeled. 

In order to include chemotaxis in our model we bias the 
spatial flux of b in Eq. (|^) by adding a chemotaxis term: 
B^hC,n{n)TiQ{9)Vn (here we take food chemotaxis as an 
example). Do(6')Vrt is the spatial derivative of the food 
concentration along the direction of movement of the bac- 
teria. Cn('^)Do(0)Vn is the derivative as sensed by the 
bacteria, where Cnl?^) can be, for example, the "receptor 
law" or a constant. The chemotaxis is attractive for 
positive values of Cn('^) and repulsive otherwise. 

The bacterial flux due to food chemotaxis translates 
in the OF model into additional flux terms in Eqs. (||- 
|). To the flux of B in Eq. (|) we add the term 
B^C,n{n) [BDil + 2 Re(pD2)] Vn (where I is the unit ma- 
trix) and to the flux of p in Eq. (^) we add the term 
B'^CnW [T>*^B+pDil]Vn . 

If a repulsive signaling material r is emitted by sporu- 
lating cells |18 l^J^ then Eqs. (§j|) are affected in a 
similar manner, with C,r{r) replacing C,n{n) and Vr re- 
placing Vn. Crif) is negative for a repulsive signal. The 
equation for the dynamics of the signaling material is 



DrV r + VrS - XbBt - XrT 



(8) 



where F^, Xb, Xr are non- negative constants. F^ is the 
rate of chemical production by sporulating bacteria, Xb 
is the rate of chemical digestion by bacteria and is the 
rate by which the chemical decompose. 

Both food chemotaxis and repulsive chemotactic sig- 
naling increase the expansion rate of the colonies. Food 
chemotaxis does not affect the colonial patterns signif- 
icantly, while chemotactic signaling does (Fig. ^). For 
short bacteria, the colonial pattern becomes less rami- 
fied, with radial branches and circular global envelope. 
For long bacteria, the global envelope becomes circular 
and the branches acquire an outward bias, changing their 
appearance from an arc-like to a hook-like appearance. 
In Fig. 1^ we also show the phenomena of global twist 
For parameters representing bacteria with inter- 
mediate length, repulsive chemotactic signaling can im- 
pose global twist on an otherwise tip-splitting pattern. 
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The twist of the branches is relative to the center of the 
colony, not to the local orientation of the branch. The 
global nature of the twist is evident by using the "de- 
chiraling" method presented by Ben-Jacob et al. in Ref. 
§ (see Fig. |d). In Ref. § Ben- Jacob et al. were able 
to obtain patterns with a global twist, using the NLD 
model and applying a rotation operator on the chemo- 
taxis term. As they argue themselves, such operator is 
inconsistent with the known biological facts and a more 
detailed model - like the OF model - is required in order 
to model any type of chirality in the bacterial colonies. 




FIG. 4. a) The effect of repulsive chemotactic signaling 
on colonies of T bacteria. To simulate short bacteria, unre- 
strained in rotation, we take Dg = 0.256 and Re(7") = 0.128. 
For chemotaxis we take (r{r) = —1, Dr = 1.0, Fr = 0.25, 
As = and Ar = 0.01. All other parameters are as in 
Due to chemotaxis the branches are radially oriented with 
circular global envelope, b) The effect of repulsive chemotac- 
tic signaling on on colonies of C. Most chemotaxis parameters 
are as in (b), Fr — 0.065 and other parameters are as in 
The branches acquire an outward bias and the global envelope 
becomes circular, c) Colony showing global twist in response 
to repulsive chemotactic signaling. Parameters values are as 
in (b), with De = 0.064, Da = 0.0625 and Re (7'') = 0.096, 
representing bacteria of intermediate length on dryer surface, 
d) "De-chiraling" the pattern in (c) by applying the mapping 
{r,<j)) + ^/R), where ii is a constant, on the polar 

coordinates (r, 0) measured from the center of the colony. 

In this letter we present a reaction-diffusion model 
which accounts for the various morphologies presents by 
colonies of P. dendritiformis. We focus on chiral features 
of the colonial patterns, but we are also able to model 
tip-splitting patterns. We successfully simulated inter- 
mediate growth patterns (not shown here). The aim of 
developing such model is to aid in the study of transition 
between the two types of growth. Our model is able to 
handle events of minute density such as a mutation in a 
single bacterium, and it will be used in following studies 
of transitions between morphologies and morphotypes. 

We thank Inna Brainis for her technical assistance. 



Presented studies are supported in part by a grant from 
the Israeli Academy of Sciences grant no. 593/95. IC 
thanks The Colton Scholarships for their support. 



[1] N. H. Mendelson, Proc. Natl. Acad. Sci. USA 75, 2478 

(1978). N. H. Mendelson and S. L. Keener, J. Bacteriol. 

151, 455 (1982). N. H. Mendelson and J. J. Thwaites, 

J. Bacteriol. 171, 1055 (1989). N. H. Mendelson, Sci. 

Progress 74, 425 (1990). 
[2] T. Matsuyama and M. Matsushita, Grit. Rev. Microbiol. 

19, 117 (1993). T. Matsuyama, R. M. Harshey, and M. 

Matsushita, Fractals 1, 302 (1993). 
[3] E. Ben- Jacob, H. Shmueli, O. Shochet, and A. Tenen- 

baum, Physica A 187, 378 (1992). E. Ben- Jacob, A. 

Tenenbaum, O. Shochet, and O. Avidan, Physica A 202, 

1 (1994). 

[4] E. Ben- Jacob et al, Phys. Rev. Lett. 75, 2899 (1995). 

[5] M. Tcherpakov, E. Ben- Jacob, and D. L. Gutnick, Int. J. 
Syst. Bacteriol. 49, 239 (1999). 

[6] E. Ben- Jacob, I. Gohen, I. Golding, and Y. Kozlovsky, 
Mathematical Models for Biological Pattern Formation, 
IMA volumes, Frontiers in Applied Mathematics Series 
(Springer- Verlag, New- York, 2000), (in press). 

[7] H. J. Gao et al, Fractals 6, 337 (1998). I. M. Sandler 
et al., Phys. Rev. E 58, 6015 (1998). I. M. Sandler et al., 
Phys. Lett. A 245, 233 (1998). 

[8] R. Weis and H. McGonnell, Nature 310,47 (1984). W. M. 
Heckl, M. Losche, D. A. Gadenhead, and H. Mohwald, 
Europ. Biophys. J. 14,11 (1986). E. M. Arnett, N. G. 
Harvey, and P. L. Rose, Acc. Ghem. Res. 22, 131 (1989). 
J. G. Heath and E. M. Arnett, J. Am. Ghem. Soc. 114, 
4500 (1992). XiaoMin Yang et al, Phys. Lett. A 193, 
195 (1994). D. VoUhardt, G. Emrich, T. Gutberlet, and 
J.-H. Fuhrhop, Langmuir 12, 5659 (1996). R. Kam and 
H. Levine, Phys. Rev. E 54, 2797 (1996). 

[9] J. M. D. Goey, G. Hinds, and M. E. G. Lyons, Europhys. 
Lett. 47, 267 (1999). 
[10] E. Ben- Jacob, Gontemp. Phys. 38, 205 (1997). E. Ben- 
Jacob, I. Gohen, and H. Levine, Adv. Phys. 49, 395 
(2000). 

[11] I. Golding, Y. Kozlovsky, I. Gohen, and E. Ben-Jacob, 

Physica A 260, 510 (1998). 
[12] I. Golding, I. Gohen, and E. Ben- Jacob, Europhys. Lett. 

48, 587 (1999). 
[13] E. Ben- Jacob et al., Physica A 282, 247 (2000). 
[14] S. Kitsunezaki, J. Phys. Soc. Jpn 66, 1544 (1997). 
[15] I. Gohen, M.Sc. thesis, Tel-Aviv University, Israel (1997). 
[16] Y. Kozlovsky, I. Gohen, I. Golding, and E. Ben-Jacob, 

Phys. Rev. E 59, 7025 (1999). 
[17] I. Gohen et al.. Fractals 7, 235 (1999). 
[18] E. Ben- Jacob et al. Nature 368, 46 (1994). 
[19] I. Gohen, A. Gzirok, and E. Ben-Jacob, Physica A 233, 

678 (1996). 

[20] J. D. Murray, Mathematical Biology (Springer- Verlag, 
Berlin, 1989). 



4 



